Appendix. Properties of multiplier matrices 

In the stability analysis of random dynamical processes we have used the properties 
of multiplier matrices following from the conditions 1 and 2 of sec. 2. Here we shall 
consider some additional properties of these matrices related to other conditions. For 
simplicity, we put to = dealing with M{^,t) = M{^,t,0) . Then eq.(9) takes the form 

-M{U) = L{y,U)M{U), 

M(e,0) = l. (A.l) 

Suppose that the eigenvalue problem for the multiplier matrix has a solution with the 
eigenfunctions forming a complete orthonormal basis: 



C/3 



g: ^i{^,t)ip^{u) = 6mn, (A.2) 



$:^„(e,%+(^,t) = i. 

n 

Property 1. If a multiplier matrix satisfying conditions (A.l) and (A.2) is Hermitian, 
then its eigenvalues are 

/i„(^,t) =exp|^ Ln{y,i,t')dt' 



O 

■ where 

? : L^{y. t) = ^+(e, t)L{y, t)^n{^, t). 

O ■ Proof. Differentiating the normalization equality in (A.2), we get 

I Using this, we find 



U / , \ , UIVI 



Since, by assumption, the multiplier matrix is Hermitian, = M , its eigenvalues are 
real, iJi*^ = Hn ■ Because of this, we have the equality 



dt "^'^ dt 
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The latter, together with (A.l) and (A. 2), gives 

from where we obtain Property 1. 

Note that L„ here is not necessarily an eigenvalue of the Lyapunov matrix, but just 
a matrix element of the latter with respect to the eigenfunctions of the multiplier matrix. 
These eigenfunctions of M are not the eigenfunctions of L . 

Another condition which yields an explicit relation between the multiplier and Lya- 
punov matrices is the commutation relation 

rt 



0, (A.3) 



Liy,U), I Liy,U')dt' 
Jo 

which sometimes called the Lappo-Danilevsky condition [26]. 

Property 2. Let the multiplier matrix enjoy conditions (A.l), (A. 2), (A.3) and be 
Hermitian, then it possesses a common set of eigenfunctions with the Lyapunov matrix, 
their eigenvalues being related as 

IIn{i,t)=expy\n{^,t')dt'Y 

where A„ is an eigenvalue of the Lyapunov matrix, 

L{y, ^, t)iPn{^, t) = A„(^, t)ipn{^, t). 

Proof. Prom (A.l) and (A.3) it follows that 

M{^,t)=expy^L{y,^,t')dt'Y 

so that [M, L] — . Commuting matrices have a common set of eigenfunctions. Because 

of the self-adjointness of the multiplier matrix. Property 1 holds true. Since now are 
the eigenfunctions of the Lyapunov matrix, we have L„ = A„ , which accomplishes the 
proof. 

The properties considered can simplify the consideration. However, a principal diffi- 
culty remains, related to the fact that the eigenfunctions ipn{i,t) depend, generally, on 
the random variable ^ as well as on time t . 
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Property 3. Assume that the multipher matrix satisfies condition (A. 2), then the 
variation (11) yields 

where 

CniC,t)^viiC,t)yiC,t), 
Proof. The proof is straightforward. 

The latter property shows that, because of the dependence of fn{^, t) on the random 
variable, one cannot, in general, define the simple relation (26) with the averaged local 
multiplier (25). This multiplier can be defined if the eigenfunctions (pn{^, t) are stochas- 
tically invariant, that is (finiC^ — fnit) , or if random fields are weak and the stochastic 
invariance of (fin{t) is an acceptable approximation. In such a case the equation from 
Property 3 reduces to 

where Cn{t) and Unit) are the same as in (23) and (25) and fn{t) — '^n{t)f . 

However, in the obtained variation the multiplier iXn{t) does not play the role of a 
quantity characterizing stability with respect to the variation of initial conditions, since 
fn{t) depends on time. Only if the eigenfunctions (pn{t) are stationary, i.e. (fn{'t)—'Pn, 
then we return to (26), and the local multiplier (25) becomes a genuine characteristic of 
local stabihty. This is why the condition of stationarity (condition 1 of sec. 2) has been 
essential for the stabihty analysis. 
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Abstract 

A particular type of random dynamical processes is considered, in which the 
stochasticity is introduced through randomly fluctuating parameters. A method 
of local multipliers is developed for treating the local stability of such dynamical 
processes corresponding to infinite-dimensional dynamical systems. The method 
is illustrated by several examples, by the random diffusion equation, random wave 
equation, and random Schrodinger equation. The evolution equation for the density 
matrix of a quasiopen statistical system subjest to the action of random surrounding 
is considered. The stationary solutions to this equation are found to be unstable 
against arbitrary small finite random perturbations. The notion of random struc- 
tural stability is introduced. 
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1 Introduction 



The stability of dynamical systems described by ordinary differential equations can be 
characterized by Lyapunov exponents. There are several kinds of the latter, for example, 
the classical Lyapunov exponents characterizing asymptotic stability [1,2], the finite- 
interval exponents describing stability on finite intervals [3,4], the pointwise local-time 
Lyapunov exponents giving information on the local stretching and contraction rates of 
trajectories [5,6], and generalized Lyapunov exponents taking into account correlations 
[7]. 

More difficult is the situation with the so-called infinite-dimensional dynamical sys- 
tems modelled by partial differential equations. There is no general stability theory for 
such equations, although in many cases one can analyse the asymptotic stability by lin- 
earizing these equations about a stationary solution [8] . 

Even more complicated is the case when partial differential equations contain ran- 
dom fields. The present paper addresses just this type of equations describing random 
dynamical processes. More precisely, only one particular type of these processes will 
be considered here, when the stochasticity is introduced into differential equations by 
means of randomly distributed parameters. Such a kind of equations may be called para- 
metrically random differential equations. This is a specific case of stochastic differential 
equations. The latter are usually treated as containing Gaussian stochastic fields related 
to the Wiener or Ornstein-Uhlenbeck processes [9-12]. These will not be touched on in 
the paper - we shall deal only with parametrically random equations. And, although the 
words random and stochastic mean the same, we shall rather use the former instead of 
the latter in order to distinguish our case from what one usually implies under stochastic 
processes. 

Parametrically random equations form an important class of equations that can be 
often met in applications, such as quantum electronics. Thus, the infiuence of interference 
on propagating signals may be modelled by randomly fiuctuating amplitude and phase. 
The random shift of frequency describes various local defects leading to the so-called 
nonuniform broadening. Some local interactions can be treated as random forces whose 
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existence may become incomparably more essential for a nonequilibrium system than the 
presence of the Nyquist noise. For example, such random local fields, due to nonsecular 
dipolc interactions, trigger pure spin superradiance in nonequilibrium nuclear magnets 
[13.14]. 

The approach suggested below may be generalized in several aspects. However, for 
the sake of clarity, here I prefer to stick myself to the parametrically random dynamical 
processes. To treat the stability of the latter, the method of multiphers is developed. For 
ordinary differential equations without random fields the method of multipliers is com- 
pletely equivalent to that of the Lyapunov exponents. However, for partial differential 
equations, especially with random fields, the method of multipliers is more natural and 
general: one may define multipliers when it is difficult or even impossible to introduce 
Lyapunov exponents. Multipliers characterize both the local as well as asymptotic stabil- 
ity. They can describe not only the exponential stability, as the Lyapunov exponents do, 
but also other types of stability. The effectiveness of using multipliers for analysing the 
local stability of dynamical processes has been demonstrated earlier for approximation 
cascades and approximation fiows [15-17]. 



2 Method of local multipliers 

Let X E B) be a set of continuous variables pertaining to a domain D which can be 
bounded or not; let t e R+ denote time; and let ^ be a set of random variables with a 
probability distribution p(^) . Consider a set, enumerated by the index i — 1,2, . . . ,d , 
of real or complex functions yi{x, ^, t) of these variables. Define the column 

y{tt)^[yi{x,C,t)] 

with respect to i and x . This column is assumed to be a solution of the evolution 
equation 

^^y{U)^v{y,U) (1) 

with the velocity field 

v{y,^,t) = [vi{x,y,^,t)] 
4 



which may contain differential and integral operators acting on y{C,t) . The evolution 
equation (1) is to be supplemented by the initial condition 

y(e,o) = / = [/,(x)], (2) 

in which fi{x) are given functions, and by boundary conditions. The quantity of final 
interest is the solution 

y{t)^ Jy{^,t)p{OdC (3) 

averaged over random variables. 

Equation (1) is a compact form of writing down a wide class of evolution equations. In 
particular, it includes the finite d -dimensional dynamical systems described by ordinary 
differential equations. But in general, this class embraces infinite-dimensional dynamical 
systems with various partial differential equations, such as can be met in considering 
hydrodynamic and plasma turbulence [18,19], dynamics of structures in liquid flows [20], 
soliton dynamics in condensed matter [21], electromagnetic radiation in large systems 
[22-24], and phase-ordering kinetics in quenched systems [25]. 

The stability of motion is related to the variation of the solution y{^,t) at time t , 

Sy{C,t) = M{^,t,to)Sy{tto), (4) 
with respect to its variation at time to . Here the multiplier matrix 

M{C,t,to)^[Mij{x,x',C,t,to)] 
is introduced with the elements 

Prom this definition it follows that 

Mij{x, x', ^, t, t) — SijS{x — x'). 

Therefore, at the coinciding times, to — t , the multiplier matrix is the unity matrix 

M{^,t,t)^i = [5ij5{x-x')] (6) 
5 



in the space of the indices i and variables x . Remind that in (4) the matrix notation 
is used according to which the right-hand side of (4) is the column 



M{C,t,to)Sy{U) 



k 

Prom the variational-derivative property 



J2 j Mik{x,x',^,t,to)5yk{x',^,to)dx' 



Syj{x',^,to) 6yk{xi,^,ti) 6yj{x',^,to) 

we find that 

M(e, t, h)M{C, h, to) = M(e, t, to). (7) 

Putting in (7) to = t and using (6), we obtain the definition of the inverse multiplier 
matrix 

M-\^,t,to)^M{^,to,t). (8) 
The properties (6)-(8) show that the multiplier matrices form a group 

MiO = {M{C,t,to) |t,toeM+}, 

which may be called the multiplier group. It is also evident that the multiplier matrices 
are the evolution operators for the variation Sy{^, t) . 

If the random variable ^ were fixed, then the transformation (4) would be contracting 
provided that ||M(^,t, io)|| < 1 ■ However, this does not yield, in general, the stability 
of the averaged solution (3). 

Variating the evolution equation (1), we get the equation 

^M(e, t, to) = e, t)M{i, t, to) (9) 
for the multiplier matrix, where the matrix 

L{y,i,t) = [Lij{x,x\y,i,t)] 

consists of elements 

U^{x,x ,y,i,t)= 10 
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In particular cases, the real parts of the eigenvalues of the matrix L{y, ^, t) define the 
Lyapunov exponents, because of which we shall call it the Lyapunov matrix. The property 

(6) plays the role of the initial condition for eq.(9). The boundary conditions for eq. (9) 
are to be obtained by variating those of eq.(l). For example, if the functions yi{x,^,t) 
are given at the boundary of D , denoted by 9 D , then 

Mij{x,x',i,t,to) = (xedB), 

and if dyi/dx is given at 9 D, then 

-^Mij{x,x',^,t,to) = Q (xedB). 

Generally, averaging the variation (4) over the random variables does not lead a simple 

relation between the variation of (3) and an averaged multiplier matrix, except for the 
case to = , when y{^, 0) is fixed by the initial condition (2) not containing ^ . In such 
a case the variation of (3) is 

5y(t) = M{t)5f (11) 

with the averaged multiplier matrix 

M{t) = I M{U,0)p{Od^. (12) 

The transformation (11) is contracting provided that | |M(t) 1 1 < 1 , from where | |52/(t) 1 1 < 
\\5f\ \ . Then we may say that the motion is locally stable on the time interval [0,t] . 

As we see, to characterize the stability of motion we need to know the multiplier 
matrix. If we are interested only in sufficient conditions of stability, we can evaluate 
the norm ||M(t)|| by employing various inequalities [26] valid for the solutions of the 
matrix equations of the type (9). Since the multiplier matrix satisfies eq.(9) involving 
the Lyapunov matrix, the properties of these matrices are closely interrelated. Some 
properties of the multiplier matrices are considered in Appendix. Throughout the paper 
we shall mainly deal with the case when the Lyapunov matrix meets the following two 
conditions greatly simplifying the analysis. 

Condition 1. The Lyapunov matrix is stationary: 

L{y,U)-m = [Lij{x,x',0]: (13) 
7 



that is, does not depend on time. 

Condition 2. The Lyapunov matrix possesses a set of eigenfunctions forming a 
complete stochastically invariant basis: 

LiOV^n = KiO^n: = [V^ni{x)], (14) 

which means that the eigenvalue problem (14) has solutions with the eigenfunctions form- 
ing a complete basis independent of the random variable ^ . 

If these conditions are not fulfilled for eq.(9), it may happen, nevertheless, that there 
exists a transformation reducing (9) to an equivalent equation with an effective Lyapunov 
matrix enjoying conditions (13) and (14). Another important case is when these conditions 
are fulfilled approximately, thus, defining an initial approximation for perturbation theory. 

Assuming condition 1, one immediately gets the solution 

M(e,t,to)=exp{L(0(t-to)} (15) 

to the matrix equation (9). Prom here, under condition (2), it follows that the multiplier 
and Lyapunov matrices possess a common set of eigenfunctions, 

M(^, t, to)(Pn = A*n(C, t, to)(pn, (16) 

with their eigenvalues related as 

l^niC, t, to) = exp{A„(0(i - to)}. (17) 

The eigenvalues of the multiplier matrix will be called the local multipliers. As is obvious 
from (17), 

fln{^,t,t)^l. (18) 

With a complete basis {(fin} , we may define the expansions for the solution 

y{Lt)=J2cn{i^t)^n (19) 
n 

and for its initial value 

f = T.fnVn, (20) 
n 
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where, according to the initial condition (2), 

The coefficient functions Cn{^,t) in the expansion (19) have the meaning of natural 
variables, or natural components, of the solution y{^,t) . For these components, from 
eq.(4) we find the variation 

SCn{C,t) = fln{^,t,to)SCn{C,to). (21) 

The expansion for the averaged solution (3) reads 

y{t)^J2<'n{t)Vn (22) 
n 

with the averaged components 

Cn{t)= J Cn{^,t)p{OdC. (23) 

Averaging eq.(16), we come to the eigenproblem 

M{t)ipn = l^nit)^n (24) 

showing that the averaged multiplier matrix (12) has the eigenvalues 

/x„(t)= J ^in{U,0)p{Odt (25) 

According to (18), we have /x„(0) = 1 . 

Using (22)-(25), we come to the conclusion that the variation (11) is equivalent to the 
set of variations 

SCn{t) = fln{t)Sfn (26) 

for the natural components (23). The apparent form of (26), in which the averaged 
multiplier (25) is just a function, makes it possible to classify the local stability properties 
at time t . 

We shall say that the n -component at time t is locally stable, locally neutral, or 
locally unstable against the variation of initial conditions if, respectively, 

\lin{t)\ < 1 {locally stable), 
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|A*n(^)| = 1 {locally neutral) , (27) 

> 1 {locally unstable). 

The motion as a whole, or the process, will be called locally stable, locally neutral or 
locally unstable according to whether sup„ |A*n(^) | is less than one, equal to or more than 
one. 

If the property of local stability, neutrality or instability of an n -component holds at 
each point t of an interval [^1,^2] , then we say that the n -component is, respectively, 
uniformly stable, uniformly neutral or uniformly unstable on this interval. Similarly, the 
process is uniformly stable, uniformly neutral or uniformly unstable on an interval [ti, t2\ , 
if the corresponding property for sup„|//„(t)| holds at each point t of that interval. 

In the case of local neutrality of an n -component the latter may be found to be unsta- 
ble with respect to higher-order multipliers defined as second, third or higher variational 
derivatives. 

As time i — > 00 , the n -component can be asymptotically stable, neutral or unstable 
if, accordingly, 

A*n(^) {asymptotically stable), 

e (0, 00) {asymptotically neutral) , (28) 

\f^n{'t) \ —>■ 00 {asymptotically unstable). 

The process is asymptotically stable, asymptotically neutral or asymptotically unstable 
if sup„ |A*n(^)| tends to zero, remains finite or tends to infinity as i — > 00 . Note that 
limt_>oo |A*n(^)| may not exist in the case of asymptotic neutrality. 

Asymptotic stability can be of different types depending on the law by which /i„(t) 
tends to zero as t ^ 00 . For instance, this can be exponential decrease, power-law decay 
or another tendency. Recall the possibility of the polynomial asymptotic stability [27] of 
solutions to stochastic differential equations. Therefore, the multipliers allow a more 
refined description of stability than the Lyapunov exponents, since the latter describe 
only the exponential asymptotic stability. 

Moreover, in the case of random processes the averaged multiplier (25) is, actually, 
the sole natural characteristic for defining stability. This follows from the fact that the 
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stability properties (27) and (28) are directly related to the multiplier fj,n{t) , but not 
to the eigenvalues of the Lyapunov matrix A„(^) , these two quantities being connected 
through the integral 

finit) = J exp{A„(OOp(Orfe (29) 

One can, of course, introduce effective Lyapunov exponents, such as effective finite- 
interval exponents 

I — to 

where 

effective pointwise exponents 

A^^^(t,t)= hm i^ln|/x,(t + At,t)| 
and effective asymptotic exponents 

A^^/= hm X^Jf{tM- 

t— »oo 

However, these would be excessive unnecessary definitions. 

What useful we can get from defining the effective Lyapunov exponents is the following. 
Assume that {A„} is a set of Lyapunov exponents ordered in the nonincreasing law, that 
is A„ > A„+i . Let be the largest integer for which I]^=i An > . Then the 
Lyapunov dimension [28] is 

1 ^ 
Dl = N + - 1 5^A„. 

I'^iV+ll n=l 

If we denote by A+ positive Lyapunov exponents, then the metric entropy is /i = A^ . 
These definitions can be transformed to the language of multipliers as follows. 
Let the local multiphers (29), at each fixed t , be ordered so that 

|A*n(^)| > |Ain+l(^)|- 

And let N — N{t) be the largest integer for which 

m 

n Mt)\>i. 

n=l 
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Then we define the local Lyapunov dimension as 



The asymptotic Lyapunov dimension is 

Dl = hm DL{t), 



t— >oo 



if this hmit exists. Similarly, the local metric entropy is 

h{t) = ]lnll\^iUt)\ {Uit)\>l)- 

^ n 

Concluding this section, it is worth making several remarks. First, the method of local 
multipliers seems to be a natural and convenient tool for analysing the local stability of 
motion. The appearance of local instability, in the case of quantum-mechanical models, 
may be connected with the quantum chaos [29] which is a temporal phenomenon. Note, 
however, that there are models of quantum systems [30,31] exhibiting the same asymptotic 
chaos as that in the models of classical mechanics. 

Second, instead of emphasizing the time dependence, one could separate out one of 
the space coordinate considering the stability properties along the chosen space direction 
[32,33]. Then we could introduce the local, with respect to the separated variable, mul- 
tipliers. Also, one could consider the motion treating a coupling parameter as a variable 
[34,35]. Then the stability properties with respect to the change of the couphng parameter 
could be analysed. 

Third, we could be interested in the stability of stochastic processes against the vari- 
ation, not of initial conditions, but of boundary conditions or of model parameters. Then 
a similar scheme of analysing the local stability could be developed. 



3 Examples of random processes 

The method described in the previous section will be illustrated here be several simple 
examples of random processes. We shall consider two types of the probability distribution 
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p(^) for the stochastic variable ^ : the uniform distribution 

p(0 = ^[e»(e + 6)-e(e-eo)], 

in which 6(^) is the unit-step function; and the Gaussian distribution 
The average value of ^ in both cases is zero, and the average of , respectively, is 

-oo o J—oo 

3.1 Random diffusion equation 

Consider the diffusion equation 

in which x E [0,L], t >0 , the diffusion coefficient D > , and ^ describes random 
fluctuations of the diffusion coefficient. Such a situation can occur, for example, in a 
diffusion process through a nonhomogeneous medium consisting of randomly distributed 
regions with different diffusion coefficients. Eq.(31) is supplemented by an initial condition 

and boundary conditions 

y(0,^,i) = ci, y{L,^,t)^C2. 
Comparing (1) with (31), we see that the velocity fleld is 

v{x, y, e, t) = iD + 0^y{.x, e, t). (32) 
For the Lyapunov matrix (10) we get 

L{x,x',0-{D + 0^5{x-x'). (33) 
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The equation (9) for the multipher matrix is 



where M — M{x, x', ^, t) . Variating initial conditions and boundary conditions gives 

M{x,x',^,Q) = d{x-x'), 

M{Q,x',^,t)=M{L,x',^,t) = 0. 

Since the Lyapunov matrix (33) is stationary, it has common eigenfunctions with the 
multipher matrix. These eigenfunctions, satisfying the boundary conditions 

iPniO) = ML) = 

and the normalization 

pL 

Jo 

are 

where n = 1,2, ... . The basis {ipn{x)} is stochastically invariant, that is, does not 
depend on the random variable ^ . 
The eigenvalues of (33) are 

A„(0 = -(i? + 0^n- (35) 

Therefore, the local multipliers (17), with to = , become 

fin{U) = exp{-{D + Oklt}. (36) 

The same answer can be obtained by solving the equation (34) yielding the multiplier 
matrix 

oo 

M{x, x', e, i) = E /""(e, t)ipn{x)ipn{x'). (37) 

n=l 

The eigenvalues of (37) are, evidently, given by (36). 
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The multiplier matrix (37) can also be found from the direct variation of the solution 

oo 



n=l 



to eq.(31), where 



L Jo 

r \ C2 - ci 
yco{x) = Ci H X. 

Now, let us find the averaged multiplier (25). For the case of the uniform noise, the 
latter is 

Mt) = ^^^^^axp(-Di;t). (38) 
The asymptotic behaviour of (38) at small and large times is 

fin{t) ~ 1 - Dklt {t ^ 0), 

^nit) ^ ^)^^^> (,^ 00). (39) 

According to (27) and (28), the behaviour of (38) defines the local stability at time t 
of an n -component. Thus, the n -component for < D is uniformly stable for all 
t > . As i — > oo , it is exponentially stable. Since sup„ |//n(i)| < 1 for all n > 1 , the 
process is uniformly stable for all t > . 

When ^0 = D , the n -component is also uniformly stable for all t > . But as 
i ^ oo , it displays now, not exponential, but power-law stability, 

/i„(t) ~ {2^oklt)-^ (t ^ oo). 

The process is uniformly stable for t > . 

If > D , then the n -component is locally stable on the interval (0, tn) , where 
tn is given by the condition \fin{tn) \ = 1 • After this, it becomes unstable for all t > tn ■ 
Because of the limit lim„^oo |A*n(^)| — oo , the process is uniformly unstable for any 
t>0 . 

In the case of the noise with Gaussian distribution, the averaged multiplier (25) be- 
comes 

finit) = exp {^kltf - Dklt^ . (40) 
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In this case, the n -component for any 7 > is locally stable till i„ = 2D/{'^knY 1 
when |//„(in)| — 1 and it is uniformly unstable for t > tn ■ As far as sup„ |A*n(^)| — 00 , 
the process is uniformly unstable for all t > . 

3.2 Random wave equation 

In the wave equation 

where x e [0, L], i > , the inclusion of the random variable ^ models fluctuations of 
the sound, or light, velocity c in a randomly inhomogeneous medium. Write the initial 
conditions to (41), 

where Ut means the derivative with respect to t . Take the boundary conditions as 

K(0,e,i) = 0, «(L,e,t) = 0. 
To reduce (41) to the standard form (1), define 

du 

Then (41) is equivalent to the system 

The corresponding initial and boundary conditions are 

yi{x,^,0) = fi{x), 1/2(2;, {,0) = f2{x), 

Z/i(0,e,i) = 0, 7/i(L,e,t) = 0. 
For the velocity field we have 

vi{x,y,^,t) = y2{x,^,t), 
V2{x,y,^,t) = {c + if^yi{x,i,t). (43) 
16 



Thence, for the Lyapunov matrix (10) we get 



L{x,x',^) 







\2 92 



Variating (41), we obtain the equations 



S{x — x') 



(44) 



9Mi 



11 



dt 

dMi2 



= M2l, 



(9M 
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M, 



dt 

dM2i 



ic + 



12 



dx'^ 



dt dt ' "^^ 

for the multipher matrix Mij — Mij{x, x',^,t) with the initial and boundary conditions 

Mij{x, x', ^, 0) = SijS{x - x'), 

Mij{0, x', e, t) = 0, Mij{L, x', e, t) = 0. 

Eigenfunctions of the Lyapunov matrix (44), satisfying the boundary and normahza- 
tion conditions 

<^m(0) = 0, iPni{L) = 0, 

["^ \iPmix)\^dx ^ 1 (i = l,2) 
Jo 

are 



(46) 



I z sm 



% 

UJr, 



where 



Tin 



UJn = (^niO = + O^n, K = -T (n=l,2,...). 



Note that these eigenfunctions are not orthogonal to each other, since 



(jj^ — ]_ 



The eigenvalues of (44) are 



(46) 
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Due to the stationarity of the Lyapunov matrix (44) , the multipher matrix has the same 
eigenfunctions with the eigenvalues 



(47) 



Equations in (45) can be solved explicitly giving 



M{x, x' , ^, i) = 2 ^ s\n{knx) sin(/c„a;') 



n=l 



cos u„t 



sin LkJrt t 



—UJn sin UJnt COS LOnt 



(48) 



It can be checked that (48) has, really, the eigenfunctions (fni{x) and the eigenvalues 
(47). The multipher matrix (48) could be also obtained from the direct variation, with 
respect to initial functions fi{x) and f2{x) , of the solution 

oo 

u{x, ^, i) = ^ [An cosa;„(^)i + sina;„(^)i] sin knX, 



n=l 



in which 

2 

= T / /i(^) sm{knx)dx, 
L Jo 

2 

^ri = TTTT / f2{x)sm{knX)dx. 

UJn{^)L Jo 

In the considered case, the eigenfunctions ipniix) depend on C, through Ci;„(^) . 
Therefore, eq.(26) is not valid. However, if the noise is weak, that is if <^ c or 7 ^ c , 
then, we can resort to perturbation theory with the zero-order basis {(fni{x)} in which 

— <^n(0) . Then the spectrum cUniC) is the first-order approximation. Respectively, 
the averaged multiplier (25) can be defined as a first-order approximation. 

Averaging (47) with the uniform distribution, we get 

l^ni{t) = ^^"fy."^^ exp{icknt), 

liol^nt 

IXn2{t)^lxl,{t). (49) 

Since |/x„i(t)| < 1 for all n > 1 and t > , the process is uniformly stable, its 
asymptotic stability being of power-law type. 
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For the Gaussian noise we find 



Hni{t) = exp ^ickj - ^ i^kjf^ , 



/^n2W=<lW. (50) 

The process is uniformly stable for all t > , and as i ^ oo , it is exponentially stable. 

Note that perturbation theory can be employed for defining the spectrum of the Lya- 
punov matrix. But after averaging over stochastic fields, the local multipliers (49) or (50) 
cannot be expanded in powers of or 7 because now these parameters enter being 
factored by the time t . 

3.3 Random Schrodinger equation 

Take the time-dependent Schrodinger equation 

for a complex wave function ip — il){x,^,t) , where x e (— oo,+oo), t > 0; m and 
U are real constants, and the random variable ^ imitates random fluctuations of the 
potential U . An initial conditions to (51) is 

A special case of the initial condition will be considered here, when the initial function 
is periodic, 

fix + L)=f{x). (52) 

We opt for a periodic initial condition in order to compare the influence of random fields 
on the Schrodinger equation with the influence of nonlinearity. When periodic boundary 
conditions are enforced for the cubic Schrodinger equation, then its periodic solutions 
arc well known to be subject to modulational long-wavelength instability [36,37]. The 
Cauchy problem for the linear Schrodinger equation with periodic initial data, but without 
random fields, has also been intensively studied [38,39]. The behaviour of the solution. 
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being closely related to incomplete Gaussian sums [39,40], has been found to look so 
chaotically as it would display the property of quantum chaos [41]. 
In the case of eq.(51), for the Lyapunov matrix we have 



S{x-x'). 



2m dx"^ 

Its eigenf unctions, satisfying the periodicity and normalization conditions 



rL 

iPn{x + L) = ipn{x), / \ipn{x)fdx = 1, 

Jo 



are the plane waves 



zfCfiX 



where n is any integer. The eigenvalues of (53) are 

The Lyapunov matrix (53) is stationary, therefore the local multiplier is 

/x„(C, t) = + + t| . 

The multiplier matrix satisfies the equation 

dM .d^M 
dt dx^ 

with the initial and periodicity conditions 



-i{U + OM 



(53) 



(54) 



(55) 



(56) 



(57) 



The solution to (57) is 



M(x,x',C,0) = 6{x-x'), 
M{x + L,x',C,t)^M{x,x',C,t). 



M{x,x',^,t)= f^n{i,t)^n{x)^*n{x'), 



(5^ 



from where it is evident that (58) has eigenfunctions (54) and eigenvalues (56). The same 
matrix (58) could be obtained by a direct variation of the solution 
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to eq.(51), where 

Cn = / f{x)ipl{x)dx. 

Jo 

The eigenfunctions (54) are stochastically invariant, which permits to define the aver- 
aged local multipher (25). For the uniform noise one gets 

The process is uniformly stable for all t > with a power-law asymptotic stability. 
For the Gaussian noise, the local multiphers are 

^,r.it) = exp + [/j t - 1 (7t)'| . (60) 

So, the process is uniformly stable for any t > with the exponential asymptotic 
stability. 

If random fields are absent, then |//„(0, i)| = 1 , that is the motion is neutral. The 
solution ip{x, 0, t) is quasiperiodic with a countable number of frequencies 

+ ("^°'- 

As is known, the behaviour of many-frequency quasiperiodic solutions can be quite com- 
plicated, often reminding chaotic one. But, of course, there is no chaos here. The neutral 
process, after the inclusion of random fluctuations, becomes stable. 



4 Stability of quasiopen systems 

From the simple models of the previous section we now pass to the consideration of a 
realistic statistical system. Take the the standard Hamiltonian 

H = J ^t(f ) _ ^(r)d 7 + 

^ J ip\r)il!\'r )^(r — ~r )ip{r )ip[r)d r d r (61) 



+ 2 



of spinless particles with a chemical potential /i , a symmetric interaction potential 
$(r) = $(— r) and with field operators ■0(r ) = ■0(r , i) . Consider the density matrix 

p{7y,t)^{'ilj\7\t)ilj{7,t)), (62) 
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in which the brackets (...) mean the statistical averaging with a statistical operator 
p(0) . The diagonal element of (62) is the density of particles 

p{7,t) = p{7,7,t)>o. (63) 

The average density of particles in a system of volume V is 

^ J^p{7,t)d7^p. (64) 

In the thermodynamic limit, V oo , the right-hand side of (64) remains constant. 
Another property of the density matrix (62) following from its definition is 

p*(f,?',t)=p(f',f,t). (65) 

Differentiating (62) with respect to time and invoking the Heisenberg equations for the 
field operators, one gets the evolution equation relating (62) with the two-particle density 
matrix 

P2(r 1, r2, r[, t) = (V'^(r'jV(^2)V'(^2)V'(n)), (66) 

in which all field operators contain the same time variable t , that is, i) = ip{r ^^t) . 

Suppose that the considered system is randomly open in the sense that it is subject to 
the action of random forces from surrounding whose influence can be taken into account by 
supplementing the evolution equation with a random fleld. Thus, the evolution equation 

for the density matrix (62) takes the form 



1 



p{ri,r2,^,t)+ 



2m 

+ j [$(f 1 - ra) - $( r2 - r 3)] ^2(^1, rg, r2, r 3, ^, t)d r^, (67) 

where the random variable ^ can, in general, be complex. If we put ^ = in (67), then 
we return to the usual evolution equation for an isolated system. 

Emphasize the difference between a random isolated system and a random open sys- 
tem. In the former case, one should add random flelds into the Hamiltonian (61), so such 
fields are to be Hermitian. In the latter case, random fields are to be inserted into the 
evolution equation, and they are not necessarily self-adjoint. 
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To make the evolution equation (67) closed, we need to resort to an approximation 
for the two-particle density matrix (66). Let us use the Hartree-Fock approximation 

P2(r 1, r2, r[, t) = p(r i, r[, t)p(r2, r'^, t) ± p(r i, r'^, t)p(r2, r[,t). (68) 

If the interaction potential $(r) is strongly singular, then the Hartree-Fock approx- 
imation leads to divergences. In that case, one must use the correlated Hartree-Fock 
approximation [42] , in which the interaction potential is smoothed by a correlation func- 
tion. Keeping in mind these both possibilities, we imply in what follows that $(r) is 
integrable. The evolution equation (67) becomes 



d 



I 



dt 2m 

$( r 1 - fa) - $(r 2 - r 3)] [p(r 1, ¥2, i)p(r 3, r 3, ^, t)± 
±p(^i, ^3, i, t)p(r3, r 2, C, t)] d f 3 . 



(69) 



The solution to eq.(69) has to satisfy conditions (63)-(65) and an initial condition 



p(r I, r 2,^,0) = /(ri, r2). 



(70) 



To analyse the stability of processes described by eq.(69), we need to consider the 
multiplier matrix 

Sp(ri,r2,^,t) 



M(ri, r2, r^, r^,^,t) 



Sp{ri,r2,C,t) 



(71) 



with an initial condition 



M(r 1, r 2, f 1, r 2, 0) = 5(r 1 - r J(5(r 2 - r 2) 



(72) 



For the multiplier matrix (71) we may write the evolution equation (9) with the Lyapunov 
matrix 

L(n.?..?;.?;.P.«.t)=^^%4^. (73) 



where the velocity field is 

v(fi,r2,p,i,t) 



I 

2m 



p(ri,r2,i,t)- 
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Prom (73) and (74) we find 

L(ri, r2,f'i, r2,p,C,i) = 
-iS(ri - r[)S{r2 - r'^) j 



d r- 



S(ri - ri)S{r2 - r^)- 



(74) 



—I 



^{ri - rg) - $(r2 - rg)] p(r3, r3,^,t)d - 

$(f 1 - r[) - $(r2 - f 'j] p(ri, ri,^,i)5(r2 - r'^)T 
$(fi - r2) - $(r2 - p(r2, r2,C,t)5(ri - r\). 



(75) 



Consider the case when the random field ^ is weak. More precisely, this means 
that the average of is much less than the average energy of the system. Of great 

importance is the limiting situation when the average of tends to +0 , that is, when 
the random field is arbitrary small although finite. A statistical system with the evolution 
equation (67) subject to the action of an asymptotically small random filed will be called 
randomly quasiopen. 

If in eq. (69) we put C = , then it is easy to check that the density matrix 

P(n, ^2, 0, t) = po(ri - ra), (76) 
where poif) is an arbitrary function satisfying just two conditions 

Po(^) = Po(- r), po(0)=p>0, (77) 
is a stationary solution. Then the density of particles 

p(r,r,0,t) = po(0) = p 

is uniform in real space. It is worth stressing that there is infinite number, actually, a 
functional continuum, of stationary solutions (76), for which conditions (77) are fulfilled. 

A common convention, concluded from the existence of stationary solutions (76), is 
that an isolated statistical system with any initial condition will finally, as i — > oo , tend 
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to a stationary state called the state of absolute equilibrium. The fact that such a state 
is infinitely degenerate, in the sense that there exist infinitely many density matrices with 
the same diagonal part but different nondiagonal parts, is interpreted as follows. All these 
density matrices are treated as statistically equivalent provided that they define the same 
set of statistical averages for local observables. As is obvious, the convention about the 
existence of absolute equilibrium is based on the assumption of its asymptotic stability. 

Using the method of multipliers we can check the stability of a stationary solution 
given by (76). The Lyapunov matrix with the solution (76) is stationary, i.e. according 
to the notation (13), 

does not depend on time. When there are no random fields, that is ^ = , then 
the eigenf unctions of L(0) are arbitrary functions of the type ^{r i — r2) with the 
eigenvalue A(0) = . The corresponding motion is neutral, as it should be expected for a 
mean-field approximation. In the presence of weak random fields, the eigenvalue problem 
for L(^) can be solved using perturbation theory. Then in the first-order approximation 

And for the local multiplier we have 

/x(e,t)=exp(-iet). (78) 

Consider the averaged local multiplier (25) assuming that the random variable ^ 
pertains to the complex plane. For the uniform distribution, let Re^ e [— 
/mCe [-e2,+6] ■ Then 

sin(6t) sinh(^2^) 
Mt) . (79) 

The asymptotic behaviour of (80) is 
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If ^1 > ^2 , then the motion is locally stable at small time. With increasing time the 



where /i(trec) = 0. . So, there occurs a peculiar stability echo. When = ^2 , then the 
first region of local stability is inside the interval < ^it < 3.8 ; the second, inside the 
interval 6.1 < ^it < 6.4 , and so on. There is no asymptotic stability, but the regions 
of local stability and instability change one another. Such a kind of behaviour is usually 
called intermittent. If ^1 < ^2 , then the motion at small time is unstable, but with 
increasing time it is again intermittent, when instability regions are interrupted by the 
stability echo. Here the intermittency happens with respect to time, but it may occur in 
space as well [43]. 

If the noise is Gaussian, we denote by 71 the dispersion of Re^ ; and by 72 , the 
dispersion of Im^ . Then the averaged local multiplier is 



For 7i > 72 , the motion is uniformly stable for all t > . If 71 = 72 , the motion is 
neutral. And it is uniformly unstable for any i > , if 71 < 72 . 

In this way, for both uniform as well as Gaussian noise there exist such small ran- 
dom fields that make the stationary solutions for the density matrix unstable. In the 
presence of such fields, the solution to (69) will wander between infinite many of the sta- 
tionary solutions, always remaining unstable and nonstationary. Therefore, the density 
of particles (63) will overlastingly depend on time fluctuating in real space. Accepting 
that no reahstic statistical system can be ideally isolated, but should be treated rather 
as randomly quasiopen, we come to the conclusion that such a system never reaches the 
state of absolute equilibrium. A randomlyly quasiopen system can become not more than 
quasiequilibrium, with perpetually persisting mesoscopic fluctuations [44]. 

5 Random structural stability 

An important question is how the qualitative behaviour of a dynamical system without 
random fields changes when the latter are switched on. Define by S — {Si} a set of 



stability is lost, but it is recovered around the recurrence time t. 



n (n = 1,2, ...) , 




(80) 
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parameters characterizing the distribution of random variables. For example, in the case 
of real random fields, S — for the uniform distribution, and 5 — ^ for the Gaussian 
distribution. In the case of complex random fields, 6 = {^1,^2} ^or the uniform noise, 
and S = {71,72} for the Gaussian noise. A dynamical system without random fields 
corresponds to S — , which implies that all Si — . 

Define the parametric neighborhood oi 5 — as a manifold {5i e [0, £j] \ei > 0} of 
stochastic parameters 5i pertaining to arbitrary small finite intervals [0, Si] with Si 
independent of each other. 

We shall say that the n -component is random- structurally stable with respect to a 
given random field, if there exists a parametric neighborhood of (5 = such that 

lim lim|/x„(t)| = lim lim 

t-+oo 5— >0 0— »0 c— >oo 

for any sequence of 5 from this parametric neighborhood. The limit 5^0 means 
the all 5j . The commutativity of the limits in the above equality can be expressed 
shorter as 

lim, lim = 0. 

_t— +00 5— >0J 

When this commutativity holds true for all n -components, we shall say that the dynam- 
ical process is random- structurally stable. 

Let us illustrate the notion of this stability for the examples of section 3. Thus, for 
the random diffusion equation we have 

lim lim |Ain(^)| = 

t— >oo (5— ►0 

for all n = 1, 2, . . . . And the opposite order of the limits gives 

{uniform noise) 

lim lim = < 

1^ 00 [Gausstan noise) 

again for any n > 1 . Therefore, the diffusion process is random-structurally stable with 
respect to the uniform noise and unstable with respect to the Gaussian noise. 
For the random wave equation we get 

lim \im\nni{t)\ = 1 

t— >oo 5— ►0 
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for all n = 1, 2, . . . and i = 1, 2 . While 



lim lim |/Xni(^)| — 

for both uniform and Gaussian noise, and any n > 1 and i — 1,2 . Consequently, the 
wave process is random-structurally unstable with respect to the uniform as well as to 
the Gaussian noise. 

Similarly to the previous case, for the random Schrodinger equation we find 

lim lim|/x„(i)| = 1, 

t— >oo 5— >0 

lim lim |//n(^)| = 

for all n e N and for both uniform and Gaussian noise. Thus, the process is also 
random-structurally unstable with respect to these noises. 
For the randomly quasiopen system of sec.4, we obtain 

lim lim I //(i) I = 1. 

t— >0O ^—^0 

But interchanging the limits makes the result undefined. For the uniform noise we get 
the intermittent behaviour, and for the Gaussian noise we can come to — oo, 1 or +oo 
depending on the relation between 71 and 72 . Thence, the dynamical process for a 
randomly quasiopen system is random-structurally unstable. 

The notion of random-structural stability is analogous to that of structural stability 
[28,45]. Neutral dynamical processes are known to be structurally unstable. As we have 
seen, such processes are also random-structurally unstable. 



28 



